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ABSTRACT 

We present observations of HESS J1640— 465 with the Fermz-LAT. The source is detected with 
high confidence as an emitter of high-energy gamma-rays. The spectrum lacks any evidence for 
the characteristic cutoff associated with emission from pulsars, indicating that the emission arises 
primarily from the pulsar wind nebula. Broadband modeling implies an evolved nebula with a low 
magnetic field resulting in a high 7-ray to X-ray flux ratio. The Fermi emission exceeds predictions 
of the broadband model, and has a steeper spectrum, possibly resulting from a distinct excess of low 
energy electrons similar to what is inferred for both the Vela X and Crab pulsar wind nebulae. 
Subject headings: ISM: individual (HESS J1640-465, 1 FGL J1640. 8-4634, 3EG J1639-4702, 
G338.3— 0.0) — pulsars: general — supernova remnants 



1. INTRODUCTION 

The basic structure of a pulsar wind nebula (PWN) 
is determined by the spin-down energy injected by the 
central pulsar and the interaction of the nebula with the 
interior regions of the supernova remnant (SNR) in which 
it evolves. Losses from synchrotron radiation in the neb- 
ular magnetic field, whose strength depends both on the 
nature of the injected wind and on the evolving size of 
the PWN, inverse- Compton scattering of ambient pho- 
tons by the energetic electron population within the neb- 
ula, and adiabatic expansion as the nebula sweeps up the 
surrounding supernova ejecta, all combine to determine 
emission structure and long-term evolution of the nebula. 
(See Gaensler & Slane 2006 for a review.) Multiwave- 
lcngth observations of PWNe provide crucial informa- 
tion on the underlying particle spectrum which, in turn, 
strongly constrains both the magnetic field strength and 
the stage of evolution. Of particular importance is the 
spectrum of particles injected into the PWN. While this 
is typically assumed to be a power law, it has long been 
known that the spectrum of the Crab Nebula is not fully 
consistent with such a structure; changes in the spec- 
tral index in the optical band seem to indicate inherent 
features in the injection spectrum, and the large popu- 
lation of radio-emitting particles may have a completely 
distinct origin from that of the higher energy particles. 

Complex structure in the observed spectrum of a 
PWN can originate in a number of ways that are as- 



1 Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge, MA 02138, USA; slane@cfa.harvard.edu; 
dcastro@head.cfa.harvard.edu 

2 Departamento de Fi'sica, Universidad Simon Bolivar, Valle 
de Sartenejas, Venezuela 

3 Kavli Institute for Particle Astrophysics and Cosmology, 
Stanford Linear Accelerator Center, Stanford, CA 94025, USA; 
funk@slac.stanford.edu; uchiyama@slac.stanford.edu 

4 Institut de Physique Nucleaire, IN2P3/CNRS, 15 rue 
Georges Clemenceau, 91400 Orsay, France 

5 Center for Cosmology and Particle Physics, New York Uni- 
versity, Meyer Hall of Physics, 4 Washington Place, New York, 
NY 10003, USA 

6 NSF Astronomy and Astrophysics Postdoctoral Fellow 

7 Universite de Bordeaux, Centre d'Etudes Nucleaircs Bor- 
deaux Gradignan, CNRS-IN2P3, UMR 5797, Gradignan 33175, 
France 



sociated with the long-term evolution, including syn- 
chrotron losses and interaction with the SNR reverse 
shock (Reynolds & Chevalier 1984; Gelfand, Slane, & 
Zhang 2009). In addition, recent studies of the spec- 
trum immediately downstream of the wind termination 
shock show that, at least in some cases, the injection 
spectrum itself deviates significantly from a simple power 
law (Slane et al. 2008), and particle-in-cell simulations 
of the acceleration process produce spectra that are well- 
described by a Maxwellian population with a power law 
tail (Spitkovsky 2008). Any such structure in the spec- 
trum of the injected particles imprints itself on the broad- 
band emission spectrum of the entire nebula. The result- 
ing breaks or regions of curvature in the emission spec- 
trum, and the frequencies at which they are observed, 
depend upon the energy at which features appear in the 
electron spectrum as well as the means by which the pho- 
tons are produced (e.g. synchrotron radiation or inverse- 
Compton emission). To fully understand the nature of 
the particle injection, as well as the long-term evolution 
of PWNe, it is thus crucial to study the emission struc- 
ture over the entire electromagnetic spectrum. 

Gamma-ray observations have provided an important 
window into the late-phase structure of PWNe. While 
the X-ray emission from older PWNe is relatively weak 
due to the long-term decline in magnetic field strength, 
inverse- Compton scattering of the cosmic microwave 
background (CMB), as well as other ambient photons, 
by the energetic particles in the nebula produce very 
energetic photons, providing unique discovery space for 
these systems. Even in younger nebulae, the 7-ray emis- 
sion provides a crucial probe of the particle spectrum 
in energy regions that are not accessible with other ob- 
servations. The inverse-Compton emission from pho- 
tons with a characteristic temperature T peaks at an 
energy e JC w 5j 2 kT for scattering from electrons of en- 
ergy r ym e c 2 . CMB photons scattered into the bandpass 
of the Fermi Large Area Telescope (LAT), for example, 
originate from interactions with electrons in the approx- 
imate energy range 0.25 - 4 TeV. Such particles produce 
synchrotron radiation with photon energies in the range 
e s ~ 0.03 — 7£?io eV, where £?io is the magnetic field 
strength in units of lOfiG. Depending on the magnetic 
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Figure 1. Fermi LAT image (2 - 200 GeV) of HESS J1640-465. 
The cyan circle indicates the uncertainty in the centroid of the 
Fermi LAT source, the magenta circle indicates the 95% encircled 
flux distribution of the HESS image, and the white circle indicates 
the 95% probability contour for the position of 3EG J1639-4702. 
The white contours outline radio emission from G338.3— 0.0 while 
the black contours at the center outline extended X-ray emission 
observed with XMM. A compact X-ray source detected with Chan- 
dra resides within the X-ray contours. 



field, which can range from > 100/xG or higher for young 
PWNe, down to ~ 5/iG for highly evolved systems, the 
synchrotron emission from particles in this energy range 
can be difficult to detect due to instrumentation limita- 
tions (at long radio wavelengths), absorption (in the opti- 
cal/UV range) or confusion with the bright sky Galactic 
background (in the infrared). Fermi measurements can 
thus provide a unique probe of the emission from a sig- 
nificant population of the particle spectrum in PWNe. 

HESS J1640— 465 (see Figure 1) is an extended source 
of very high energy 7-ray emission discovered with the 
High Energy Stereoscopic System (H.E.S.S.) during a 
survey of the Galactic plane (Aharonian et al. 2006). 
Centered within the radio SNR G338. 3-0.0 (Shaver & 
Goss 1970), the deconvolved TeV image of the source has 
an RMS width of 2.7 ± 0.5 arcmin (Funk et al. 2007). 
HI measurements show absorption against G338.3— 0.0 
out to velocities corresponding to the tangent point, in- 
dicating a distance of at least 8 kpc (Lemiere et al. 
2009), and thus implying a rather large size for the PWN 
(Rpwn > 6.4c?io pc, where dio is the distance in units 
of 10 kpc). X-ray observations with XMM (Funk et al. 
2007) and Chandra (Lemiere et al. 2009) establish the 
presence of an accompanying X-ray nebula and an X-ray 
point source that appears to be the associated neutron 
star. In addition, Aharonian et al. (2006) noted the pres- 
ence of the unidentified EGRET source 3EG J1639-4702 
located at a nominal position 34' from HESS J1640— 465, 
but the very large error circle on its position made any as- 
sociation with HESS J1640— 465 highly uncertain. Here 
we report on observations of HESS J1640— 465 with the 
Fermz-LAT. The observations and data analysis are sum- 
marized in Section 2, and a discussion of the observed 
7-ray emission is discussed in the context of the evo- 
lutionary state of HESS J1640-465 in Section 3. Our 
conclusions are summarized in Section 4. 
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Figure 2. Profile of the Fermi LAT emission from HESS 
J1640— 465. The histogram corresponds to the best-fit point source 
profile, including the diffuse Galactic and extragalactic background 
which is shown separately in red. 

2. OBSERVATIONS AND DATA ANALYSIS 

We investigated 7-ray events acquired from the region 
surrounding HESS J1640-465 with the Fermi-LAT dur- 
ing the period 2008 August 5 to 2009 November 13. Stan- 
dard event selection was applied, using the "diffuse" class 
as defined in Atwood et al. (2009), using events with 
zenith angles less than 105 degrees to minimize the por- 
tion of the Earth limb in the LAT field-of-view (Abdo et 
al. 2009a). The "Pass6 version 3" instrument response 
functions were used. Standard analysis tools available 
from the Fermi Science Support Center (version v9rl5p2) 
were used for the reduction and analysis. In this work, 
the mapcube file gll_iem_v02.fit is used to describe the 
Galactic 7-ray emission, and the isotropic component 
is modeled using the isotropicJem_v02.txt table. Data 
analysis details follow those in Castro & Slane (2010). 

For source detection and spatial analysis of the source, 
we used only events with energies in the range 2 - 200 
GeV converting in the front section of the tracker, where 
the Point Spread Function (PSF) is narrower. At 2 GeV, 
the 68% containment radius of the PSF is approximately 
18 arcmin, which is considerably larger than the source 
of interest. Based on an unbinned maximum likelihood 
analysis, using the gtlike routine, a LAT source coin- 
cident with HESS J1640— 465 is detected with a signifi- 
cance of ~ 17cr based on the TS statistic (Mattox et al. 
1996). The counts map is shown in Figure 1, along with 
the uncertainty in the centroid position based on our 
analysis (cyan circle). We have also included the HESS 
error circle (magenta) and contours from the XMM ob- 
servation of HESS J1640-465 (black) and the MOST 
observation of G338.3— 0.0(white contours). The error 
circle for 3EG J1639— 4702 is indicated by a dashed white 
circle^ The centroid of the LAT emission is located at 
16 h 40 m 46 s , -46°30'44", in good agreement with the po- 
sition of HESS J1640— 465, and the brightness distribu- 
tion is consistent with an unresolved source (Figure 2); 
models for a disk of extent between 0.05 — 0.2 degrees sig- 
nificantly degrade the quality of the fit, providing strong 

8 EGRET position contours are not necessarily circular; the ra- 
dius of the contour shown, from Hartman et al. (1999), contains 
the same solid angle as the formal 95% contour. 
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Figure 3. Fermi LAT spectrum of HESS J1640-465. Statisti- 
cal (systematic) uncertainties are indicated by solid (dashed) error 
bars. The dashed line corresponds to the best-fit power law model 
described in the text. 

evidence for an unresolved source. 

The source spectrum was extracted using both front 
and back events, and covering the energy range 0.2 — 
51.2 GeV. Standard background models were used to ac- 
count for both Galactic and extragalactic diffuse emis- 
sion as well as instrumental background. Contributions 
from field sources identified in the one- year Fermi-LAT 
First Source Catalog (Abdo ct al. 2010a) were in- 
cluded in the analysis. The LAT spectrum for HESS 
J1640— 465 is shown in Figure 3, and is well-described 
by a power law with T = 2.30 ± 0.09 and a F(> 
100 MeV) = (2.8 ± 0.4) x 10" 7 photons cm~ 2 s" 1 [or 

(1.7 ± 0.2)) x 10~ 10 erg cm" 2 s _1 in the 0.1 - 300 GeV 
band] based on spectral fits for which statistical error 
and systematic errors associated with uncertainties in 
the Galactic background level (estimated by artificially 
changing the Galactic background normalization by ±3% 
- see Abdo et al. 2009b) are added in quadrature. 

We note that Fermi observations of pulsars reveal spec- 
tra with distinct cutoffs above ~ 1 — 10 GeV (e.g. Abdo 
et al. 2009). The absence of such a cutoff in Figure 3 
implies that that the bulk of the emission does not arise 
directly from an unseen pulsar in HESS J1640— 465. Our 
fits imply a lower limit of 40 GeV (at the 90% confidence 
level) for any exponential cutoff energy in the spectrum. 
Joint fits with the HESS spectrum are also well-described 
by a power law. Addition of a second power law, with 
an exponential cutoff, does not improve the fit. Such a 
fit can, however, accommodate ~ 20% of the observed 
flux in the second power law for cutoff energies between 
1 — 8 GeV (beyond which such a component is statisti- 
cally disfavored). 

The source 1FGL J1640.8-4634, from the one-year 
catalog, is coincident with the source we detect, and the 
quoted flux is in agreement with our measurements as 
well. The flux of the source is consistent with that of 
3EG J1639— 4702; the spectral index is somewhat flat- 
ter, though also consistent within the uncertainties. 

3. DISCUSSION 

The evolutionary state of a composite supernova rem- 
nant system is strongly constrained by the observed size 
of the SNR and the inferred spin-down properties of the 



associated pulsar. Radio observations establish a radius 
Rsnr ~ 11.6dio pc for G338.3— 0.0. The observed extent 
of HESS J1640-465 constrains the radius of the PWN to 
Rpwn > 6.4dio pc0 For evolution in the Sedov phase, 
the SNR radius is 

iw=^0^y V s P c a) 

where E§i is the supernova explosion energy in units 
of 10 51 erg, no is the number density of the ambient 
medium, and £3 is the SNR age in kyr. The evolution 
of the PWN radius through the SNR ejecta is given by 
(Chevalier 1977) 

RpWN ~ l-l-Bo i 38^5i 10 -^i0 1 ^ 2i 3^ 5 P C ( 2 ) 

where i?o,38 is the initial spin-down power in units of 
10 38 erg s _1 and M 10 is the ejecta mass in units of 1OM . 
Figure 4 illustrates the evolution of the SNR and PWN 
radii under such assumptions for a range of values for the 
ambient density and initial spin-down power, assuming 
an ejecta mass of 8M and E51 = 1. For the observed ra- 
dius of G338.3— 0.0, we see that the age must be around 

5 — 8 kyr for reasonable ambient conditions. The PWN 
has ideally expanded to a larger radius by this time, 
which means that under real conditions the SNR reverse 
shock has already encountered and disrupted the PWN. 
This is illustrated in Figure 4 where the solid blue and 
red curves show the SNR and PWN radii (respectively) 
as a function of time for the model shown in Gelfand 
et al. (2009) with E = 10 40 erg s" 1 , M ej = 8M , 
no = 0.1 cm -3 , and E 5 i = 1. Here the early SNR 
evolution is calculated using the solution from Truelove 

6 McKee (1999), and the SNR and PWN evolution is 
treated self-consistently. The SNR curve initially rises 
more quickly than the dashed blue curves, which assume 
a Sedov solution from the onset, but approaches the Se- 
dov solution at later times. The solid red curve shows a 
distinct reduction in the radius of the PWN upon com- 
pression by the SNR reverse shock. Such a reverse shock 
interaction is consistent with both the inferred size of 
HESS J1640— 465 and with the observed offset between 
the putative pulsar and the surrounding nebula (Funk et 
al. 2007, Lemiere et al. 2009). 

A full exploration of the parameter space constrained 
by both the spatial and spectral properties of HESS 
J1640-465 and G338.3-0.0 is beyond the scope of this 
initial investigation. Here we have modeled the PWN 
emission following the description in Lemiere et al. 
(2009). A simple power law injection of particles from 
the pulsar into a 1-zone nebula is assumed, and radiative 
losses are assumed to be dominated by synchrotron ra- 
diationJ3 The time evolution of the pulsar spin-down is 

9 Here, and throughout, Rsnr refers to the distance from the 
SNR center to the outer edge of its blast wave, while Rpsr refers 
to the distance from the pulsar to the outer boundary of its wind 
nebula. 

10 This assumption is valid for magnetic fields larger than ~ 
1/xG, which is typically the case in PWNe. We note, however, 
that at late times (ages of order 10 kyr, as suggested for HESS 
J1640— 465), when the magnetic field decreases due to expansion, 
IC losses can begin to be significant. 
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Figure 4. Time evolution of the SNR and PWN radii for a range 
of values for the ambient density and initial spin-down power of the 
pulsar. The solid curves correspond to models from Gelfand et al. 
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where tq is a characteristic spin-down timescale for the 
pulsar and b is the pulsar braking index. Pulsations from 
the putative pulsar have not yet been detected. To esti- 
mate the current spindown power, we use the empirical 
relationship for E/L x obtained by Possenti et al. (2002), 
which yields E = 4x 10 36 erg s _1 . Based on the observed 
scatter in the E/L x , the uncertainty on this estimated 
value may be a factor of 10 or larger. The maximum 
particle energy is assumed to be limited by the condition 
that the particle gyroradii do not exceed the termination 
shock radius (de Jager & Harding 1992). The lower limit 
on the particle energy is a free parameter, and the overall 
normalization is set by the wind magnetization parame- 
ter a, equal to the ratio of power in particles to that in 
the Poynting flux. The magnetic field in the nebula is 
assumed to evolve as 



B{t) = B Q / [1 + (t/r ) a ] 



(4) 



where a is a free parameter. 

The broadband emission model results are shown in 
Figure 5 where we plot the Fermi and H.E.S.S spectra 
along with the radio upper limit from GMRT observa- 
tions (Giacani et al. 2008) and spectral confidence bands 
derived from Chandra (Lemiere et al. 2009). The black 
curves represent the model prediction for the synchrotron 
(left) and inverse- Compton (right) emission that best de- 
scribes the X-ray and TeV 7-ray spectra, similar to re- 
sults from Lemiere et al. (2009); the parameters for the 
model, which were adjusted "by hand" to provide good 
agreement with the radio, X-ray, and TeV 7-ray data, as 
well as the inferred size of the PWN, are summarized in 
the caption. As seen in Figure 5, this model significantly 
underpredicts the observed Fermi-LAT emission. 

As discussed above, our spectral fits can formally ac- 
commodate up to about ~ 20% of the observed flux in 
a pulsar-like component characterized by a power law 
with an exponential cutoff energy between 1 and 8 GeV. 
This corresponds to an energy flux (above 100 MeV) of 
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Figure 5. Electron spectrum (upper) and broadband emis- 
sion model (lower) for HESS J1640— 465 assuming the evolution- 
ary history described in the text. The black curves represent a 
PWN with an age T = 10 kyr, and B(T) = assuming 
E(T) = 4 X 10 36 erg s — 1 and an injection spectrum with a = 10~ 3 , 
7 = 2.5, tq = 500 yr, and E mln = 115 GeV. The magnetic field evo- 
lution is characterized by a = 0.65. The magenta curves represent 
the scenario with a low-energy Maxwellian electron component re- 
placing the low-energy portion of the electron power-law spectrum. 
The mean temperature for the IR and optical photon fields are 15 K 
and 5000 K, respectively, and the energy densities relative to the 
CMB are 4 and 1.15. The dashed curve in the upper panel repre- 
sents the truncated portion of the power law that was replaced by 
a Maxwellian. The dashed blue curve in the lower panel represents 
a model for which all of the 7-ray emission results from pion decay. 

~ 3.8 x 10~ n erg cm~ 2 s _1 . For the spin-down power 
suggested by its X-ray luminosity, the available energy 
flux from the pulsar that powers HESS J1640— 465 is 
- 3.3 x lO" 10 ^ 2 erg cm" 2 s _1 . Thus, as much as 10% 
of its spin-down could conceivably be contributing di- 
rectly to 7-rays in the LAT band. There are known 
radio pulsars within the field of 1FGL J1640.8-4634 
as well. PSR J1637— 4642, for example, is located 
within ~ 37 arcmin and has a total energy flux of 
2 x 10~ 10 erg cm~ 2 s _1 given its estimated distance and 
measured spin-down. Thus, while our spectral fits do 
not require a pulsar-like component, it it quite feasible 
that one or more of these sources contributes as much as 
20% of the observed flux. Inspection of Figure 5 makes 
it clear, however, that the remaining LAT emission still 
vastly exceeds the predicted flux from HESS J1640— 465. 

We note that simple power-law injection models for 
Vela X, another evolved PWN, fail to reproduce the ob- 
served broadband spectrum (LaMassa, Slane, & de Jager 
2009). The presence of an excess population of low- 
energy electrons has been suggested, and models for the 
inverse-Compton scattering of photons from this popu- 
lation predict an excess of 7-rays in the GeV range (de 
Jager, Slane, & LaMassa 2009). This excess has been 
confirmed with observations by both AGILE (Pellizzoni 
et al. 2010) and Fermi (Abdo et al. 2010b). Moti- 
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vated by this result, we modified the evolved power law 
spectrum from our model for HESS J1640— 465 by trun- 
cating the lower end of the power law and adding a 
distinct low-energy component. Based on results from 
simulations of shock acceleration (Spitkovsky 2008), we 
chose a Maxwellian distribution for this population. Our 
resulting (ad hoc) particle spectrum is shown in the 
upper panel in Figure 5, and the resulting broadband 
emission is shown in the magenta curves in the lower 
panel. Here we have adjusted the normalization of the 
Maxwellian to reproduce the emission in the Fermi-LAT 
band, which is produced primarily by upscattered in- 
frared (IR) photons from local dust. The energy den- 
sity and mean temperature of the IR photon field was 
adjusted slightly to improve the agreement between the 
data and the model, but the values (listed in the caption) 
are within reasonable expectations (see, e.g., Strong et 
al. 2000). We find a mean value of 7 2 x 10 5 for the 
electrons in the Maxwellian component, and roughly 5% 
of the total electron energy in the power law tail, con- 
sistent with results from particle-in-cell simulations P"1 
The associated pair multiplicity relative to the integrated 
Goldrcich- Julian injection rate is of order 10 6 , similar 
to that inferred for the Crab Nebula as well as several 
other PWNe (see Bucciantini et al. 2010). Recent work 
by Fang & Zhang (2010) uses a similar input distribu- 
tion to successfully model the emission for several PWNe 
including HESS J1640— 465. However, their results for 
HESS J1640— 465 underpredict the observed GeV emis- 
sion from this source, apparently due to use of a slightly 
lower bulk Lorentz factor and a larger fraction of the to- 
tal energy in the power law tail than we have used in this 
analysis. 

We note that the Maxwellian shape for the low-energy 
electron population is not unique. Indeed even a very 
narrow Gaussian distribution can produce the GeV 7- 
ray emission without exceeding the radio upper limit for 
the synchrotron emission. For any of these models, the 
total energy in electrons requires a larger initial spin- 
down period than assumed for the simple power-law in- 
jection models - by a factor of roughly 3 for the adopted 
Maxwellian distribution. This is within the uncertainty 
in our assumed value based on scaling from the X-ray 
luminosity of the PWN. It is also important to note that 
while our proposed model is consistent with the observed 
properties of this system, there are potential degenera- 
cies in the effects of different model parameters, meaning 
that the proposed scenario is not necessarily unique. 

An alternative scenario for the 7-ray emission is that 
it arises from the SNR itself, and not the PWN. Non- 
thermal bremsstrahlung has been suggested as a mecha- 
nism for the production of 7-rays in SNRs (e.g. Bykov 
et al. 2000). This process can be dominant for electron- 
to-proton ratios of order 1. However, the value typical 
of local cosmic rays is closer to 10~ 2 , and even smaller 
values appear favored in models for 7-ray emission from 
SNRs (e.g. Morlino et al. 2009, Zirakashvili & Aha- 
ronian 2010, Ellison et al. 2010), so that this process is 
typically not dominant. The dashed blue curve in Figure 
5 represents a model for the emission from the collision of 
protons accelerated in the SNR with ambient material, 
leading to 7-rays from the production and subsequent 

11 A. Spitkovsky, private communication. 



decay of neutral pions. The 7-ray spectrum is calcu- 
lated based on Kamae et al. (2006) using a scaling fac- 
tor of 1.85 for helium and heavy nuclei (Mori 2009), and 
we have used a power law distribution of protons with 
dNp/dE cx E~ 2A to best reproduce the observed spec- 
trum. Assuming a shock compression ratio of 4 and that 
25% of the total supernova energy appears in the form of 
relativistic protons, an ambient density no ~ 100 cm~ 3 is 
required to produce the model shown in Figure 5. This 
is much higher than can be accommodated for the ob- 
served size of the SNR and the lack of observed thermal 
X-ray emission from the SNR. Such high densities are 
found in dense molecular clouds, suggesting that the 7- 
rays could be produced by particles that stream away 
to interact with high-density material outside the SNR. 
However, only the most energetic particles can escape the 
acceleration region, which is in conflict with the proton 
spectrum we require to match the data. Moreover, the 
observed TeV emission appears to originate from within 
the SNR boundaries, making such an escaping-particle 
scenario appear problematic. Based on this, along with 
the lack of a spectral cutoff that might suggest emission 
from a central pulsar, we conclude that the GeV 7-ray 
emission most likely arises from the PWN. 

It is well-known that a distinct low-energy electron 
population resides in the Crab Nebula, although the ori- 
gin is not well-understood. Atoyan (1999) has suggested 
that the spin-down timescale r (see Equation 1) may it- 
self be time-dependent, resulting in a large energy input 
in the earliest epoch of pulsar spin-down, when signif- 
icant synchrotron and adiabatic losses would result in 
rapid cooling of the electrons. Studies of the broadband 
emission from 3C 58 indicate an injection spectrum that 
differs from a pure power law as well (Slane et al. 2008). 
As noted above, observations of Vela X appear to re- 
quire an excess of low-energy electrons which may be a 
relic population or could have been produced through 
rapid synchrotron losses associated with the increased 
magnetic field strength during the reverse-shock crush- 
ing stage. More complete modeling of the broadband 
emission of HESS J1640— 465, accounting for the full dy- 
namical evolution of the system, including the effects of 
the reverse shock on the PWN, is required to assess the 
overall energetics and underlying particle spectrum more 
completely, but is beyond the scope of the results we re- 
port here. 

4. CONCLUSIONS 

Broadband studies of HESS J1640— 465 have identi- 
fied this source as a likely PWN, with X-ray observa- 
tions providing images of an extended nebula as well as 
the putative pulsar powering the system. Modeling of 
the PWN evolution based on inferred parameters of the 
pulsar imply detectable emission in the GeV 7-ray band, 
and our investigations of the Fermi-LAT observations of 
this region reveal clear evidence of such emission, con- 
sistent with the source 1FGL J1640.8-4634. The flux 
and spectrum we derive are consistent with that of the 
previously-identified source 3EG J1639— 4702, and the 
much-improved position makes it likely that the emission 
arises from HESS J1640— 465. The lack of a spectral cut- 
off rules out an association with emission directly from 
the pulsar that powers the nebula, and the flux is in- 
consistent with 7-rays from G338.3— 0.0, in which HESS 
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J1640-465 resides. 

We have investigated the radio, X-ray, and 7-ray emis- 
sion from HESS J1640— 465 in the context of a simple 
one-zone model in which power is injected into a neb- 
ula at a time-dependent rate consistent with the current 
observed X-ray emission from the system. We find that 
models constrained by the observed size of the associated 
SNR, as well as limits on the size of the PWN, require an 
approximate age of 10 kyr and a current magnetic field 
strength of only ~ 4/xG, consistent with expectations for 
the late-phase evolution of a PWN. The conditions in 
such an evolved PWN yield a considerably higher 7-ray 
flux, relative to the X-ray flux, than in younger systems 
where the higher magnetic field results in significant syn- 
chrotron radiation. 

The observed Fermi-LAT emission from HESS 
J1640— 465 significantly exceeds that predicted by our 
broadband models. We propose that the excess emission 
is a signature of a distinct population of low-energy elec- 
trons similar to that inferred from studies of the Crab 
Nebula and Vela X, although the nature of this electron 
component is not well constrained. Deeper radio obser- 
vations are needed to place stronger constraints on this 
population. Sensitive searches for pulsations from the 
central pulsar are of particular importance to constrain 
the spin-down properties of the system, which can only 
be very roughly constrained at present. There has been 
considerable success in identifying such pulsars with the 
Fermi-hAT ', but the lack of an obvious pulsar-like spec- 
trum in HESS J1640— 465 may argue for more likely suc- 
cess with deep radio timing searches. 
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